Charging effects in a quantum wire with leads 

V. A. SablikovQ 

Institute of Radio Engineering and Electronics, Russian Academy of Sciences, Fryazino, Moscow District, 141120, Russia 



o 
o 
o 



o 

(N 
> 

in 
a^ 
m 
oo 
O 



I 

O 

o 



S. V. Polyakov 

Institute for Mathematical Modelling, Russian Academy of Sciences, Miusskaya sq.^a, Moscow, 125047, Russia 

M. Biittiker 

Departement de Physique Theorique, Universite de Geneve, CII-1211, Geneve 4, Switzerland 
(published in Phys. Rev. B, 61, No 20, 13763-13773 (2000)) 

We investigate the distribution of the electron density and the potential in a quantum wire 
coupled to reservoirs, treating this structure as a unified quantum system and taking into account the 
Coulomb interaction of electrons. The chemical potential difference that exists between a decoupled, 
isolated quantum wire and the reservoirs gives rise to charge transfer in the coupled system. We 
show that the quantum wire can be charged positively or negatively or remain neutral as a whole, 
depending on such factors as the wire radius and the background charge density in the wire. The 
magnitude of the charge and its sign are to a large extent determined by the exchange interaction 
of the electrons in the wire. Using a Hartree-Fock approach, we develop a model of a quantum 
wire which includes the reservoirs. This model allows us to find the self-consistent distribution 
of the electron density and the potential in the wire both at equilibrium and in the presence of 
transport. The linear conductance is investigated as a function of the chemical potential. The 
nonadiabatic transition from the reservoirs to the wire leads to conductance oscillations caused by 
multiple scattering of electron waves. The period of the oscillations depends on the charge acquired 
by the wire and the exchange energy. We find that the exchange interaction strongly enhances the 
Friedel oscillations near the contacts. However, they do not noticeably suppress the conductance 
because the wire has a finite length and is charged. Under far from equilibrium conditions, which 
appear when the applied voltage exceeds the Fermi energy in the wire, the system becomes unstable 
with respect to fluctuations of the electric potential and the electron density. The instability results 
in the appearance of multistable electron states. 



I. INTRODUCTION 



Electron transport in quantum wires (QWs) now at- 
tracts considerable interest because of the fundamental 
importance of the electron-electson (e-e) interaction in 
one-dimensional (ID) systems.Ela In addition, it is ex- 
pected that the e-e interaction may produce important 
effects in the transport that would be attractive for ap- 
plications. However, up to now there is no clear and 
unambiguous knowlege of which effects of e-e interac- 
tion are observable under the realistic conditions of a 
QW with leads and what the reasons are for the devi- 
ations frppi conductance quantization (observed experi- 
mentallyu), and very little is known about the electron 
transport under far from equilibrium conditions. 

In order to understand the transport properties, it is 
necessary to know the spatial distribution of the electric 
potential and the electron density in the stractures under 
investigation. Many arguments and factsQcl show that 
the leads play an essential role in the conductance, if the 
transport is investigated by measuring the electric cur- 
rent in an external circuit and by measuring the voltage 
drop between the leads. It is obvious that the interaction 
of the QW with the leads is not weak. 



In the present paper we investigate this interaction, 
considering the QW and the leads as a unified system. An 
important problem that arises is how the electric field and 
the electron density are distributed at equilibrium (i.e., in 
the ground state), and how they are redistributed when 
an external voltage is applied. We show that even in the 
equilibrium state the QW acquires a charge and a con- 
tact potential difference between the QW and the leads. 
This phenomenon is similar to the well known contact 
potential difference in classical conductor systems. The 
essential difference between the low-dimensional system 
under consideration and classical conventional 3D con- 
ductors is that the contact potential is not screened over a 
finite length but spreads over a lepigth determined by the 
geometrical size of the structure.LI The contact potential 
difference is determined by the difference in the chemical 
potentials A/i of electrons in the QW and in the leads, 
when these subsystems are considered independently. By 
analyzing the chemical potentials we find that three cases 
are possible, depending on the wire radius and the back- 
ground charge density: (i) the QW is charged positively, 
(ii) the QW acquires a negative charge, or (ui) the QW 
remains neutral as a whole. 

We develop a model of a QW with leads in which the 
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lead- wire interaction is taken into account. It is based 
on a Hartree-Fock approach for the electrons in the QW 
and the representation of leads (which are considered as 
electron reservoirs) in a way that takes into account their 
3D nature, but requires only a ID calculation. Using this 
model we investigate the distribution of the electron den- 
sity and the electric potential at equilibrium, as well as 
under far from equilibrium conditions. In particular, we 
emphasize the role of the exchange interaction effect on 
these quantities. If in the absence of coupling the chem- 
ical potential of the wire exceeds that of the reservoir, 
> 0, the coupled system expels electrons from the 
QW into the leads, when equilibrium is established. As 
a result a potential well appears in the QW. The wire acts 
as a charge donor to the reservoir. If A/i < 0, electrons 
are attracted by the wire, and a potential barrier arises 
between the QW and the leads. The wire acts in this case 
as an acceptor. In the case where Afi = 0, the electron 
density is redistributed only within the QW and only 
Friedel oscillations arise near the contacts. The Friedel 
oscillations appear in all cases; however, they are super- 
imposed on a much more slowly varying potential created 
by the acquired charge. The slowly varying potential has 
an amplitude that is large compared to the Friedel oscil- 
lations and thus dominates the scattering processes. The 
exchange interaction strongly affects the potential shape 
and somewhat enhances the Friedel oscillation amplitude. 

When an external voltage is applied, the chemical po- 
tentials in the electron reservoirs are shifted relative to 
each other, disturbing the electron flows in the QW. As 
this takes place, the electron density, the potential, and 
the exchange energy are changed self-consistently. The 
importance of the electrostatic potential distribution in 
quantum wires with leads, especially for the investigation 
of time-dependent transport, andJor nonlinear transport, 
has been emphasized previouslyJ13 but quantitative cal- 
culations have to our knowledge not been reported thus 
far. 

Under far from equilibrium conditions, when the ap- 
plied voltage exceeds the Fermi energy in the QW, the 
electron density is substantially redistributed between 
the QW and the reservoirs, giving rise to a very strong 
variation in the potential landscape. In turn the poten- 
tial produces a variation in the electron density. The 
connection between the electron density and the poten- 
tial is very important for the understanding of nonlinear 
transport. The need of a self-consistent treatment has 
been emphasized by Landauer.Eil Within the scattering 
approach, in the weakly nonlinear regime, it has been 
vestigated by Christen and oim-of the present authorsO 
and by Ma, Wang and Guo.EJ The calculation of the 
present paper allows us to investigate the strongly non- 
linear transport in the system under consideration. If 
the applied voltage is high enough, the self-consistent 
connection of charge and potential gives rise to an insta- 
bility of the electron density distribution and ultimately 
leads to multistability of the electron states in the QW. 
This means that several stable states with different spa- 



tial distributions of the electron density and the potential 
are possible at a given applied voltage. 

The paper is organized as follows. In Sec. ^ the chem- 
ical potential difference between the decoupled QW and 
the 2D electron reservoir is analyzed. Section [II de- 
scribes the model of the QW with leads. Section IV 



contains the results of the numerical calculations of the 
electron density and potential distribution in the QW 
with leads. In Sec. ^ the multistability of the electron 
states is described that appears for far from equilibrium 
conditions. 



II. CONTACT POTENTIAL DIFFERENCE 

To be specific we consider a QW connecting two re- 
gions of a 2D electron gas. We assume that there are no 
nearby gates and that all electric field lines emanate and 
terminate either on the wire or on the 2D electron gas. 
QW structures of this kind are produced by etching of 
heterostructures with a 2D electeoa-gas. Such structures 
are widely used in experiments J1j~I13 

First we investigate such a QW separately from 2D 
electron reservoirs. The uncoupled wire is charge neu- 
tral. The electron charge is concentrated inside the QW, 
while the compensating positive charge of the impurities 
is really located in the immediate vicinity of the QW or at 
its surface. The decoupled QW and the electron reservoir 
have their own chemical potentials /iiD and /i2D which are 
generally not equal each other. We are interested in the 
chemical potential difference Afi between the QW and 
the reservoir. |_. 

According to Seitz's theorenxj the chemical potential 
in the QW is determined by the Fermi energy of nonin- 
teracting electrons and the self-energy T,{kp) which takes 
into account the e-e interaction. 



Mid = eo + £_f + '^{kp) ■ 



(1) 



Here Eq « ir'^fi'^ / (2ma^) is the first subband energy 
caused by transverse confinement (a is the QW radius, 
m is the effective mass of electrons). E(fci?) contains the 
contributions arising from the exchange and correlation 
interaction as well as from the electron interaction with 
the positive background charge. The exchange and cor- 
relation energy was-|investigated in the recent paper of 
Calmels and GolclE2l using the seLLponsistent theory of 
Singwi, Tosi, Land and SjolandertiJ for the case where 
only the lowest subband is occupied. The Hartree en- 
ergy is easily estimated if we assume that the positive 
charge is located at the surface of the QW. These calcu- 
lations lead to the following expression for the chemical 
potential /iid of the QW in terms of the dimcnsionlcss 
parameters — l/{2aBn) and P = as/a (with as the 
effective Bohr radius and n the ID electron density): 
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where Ry is the efTective Rydberg. In Eq. (||) the first 
term is the lowest subband energy, the second term is 
the kinetic energy, and the third term represents the ex- 
change and correlation energy. There are two expressions 
for the exchange-correlation term depending on whether 
r, < 1 or > 1: 
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5.577r - Apr, 
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if Ts < 1 



and 



1.84 
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if r. , > 1 



The last term in Eq. (g) is the Hartree energy of the 
electron interaction with the positive background; Bh is 
a numerical factor that depends on the radial distribu- 
tion of the electron density. If the electron density is 
distributed uniformly, Bh ~ 1/3. 

The chemical potential in the 2D reservoi|t_can eas- 
ily be obtained from the known expressionE3 for the 
electron energy as a function of the density parameter 
Rs = (7ra^A^2D)~"'^ (where A^2D is the 2D electron den- 
sity). In the high density case, Rg < \/2, one obtains 
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0.0863i?Jni?, + 0.519i?, , 



(3) 



where d is the thickness of the 2D layer. 

Using Eqs. (|^) and (||), the chemical potential differ- 
ence A/z = /ziD — M2D is calculated as a function of the 
QW radius and the density parameter r,. The results of 
these calculations are illustrated in Fig. |l], where the A/i 
dependence on the wire radius a is shown for various val- 
ues of . Here we consider as an independent param- 
eter because the background charge density depends on 
external factors, such as the charge absorbed at the wire 
surface. Figure |l| shows that in sufficiently thin wires 
the chemical potential is higher than in the reservoir. 
However, with increasing radius of the QW the chemi- 
cal potential in the wire can become lower than in the 
reservoir. 

The chemical potential difference Afi is caused by all 
energy components contributing to the chemical poten- 
tials in the QW and in the reservoir. As an example it 
is instructive to consider the estimations for the specific 
case where the QW diameter is equal to the 2D layer 
thickness, d — 2a, and the background charge density 
per unit area, N2r>, is the same in the reservoir and in 
the QW. The latter means that n = 2aA'2D- In this 
case the energies contributing to the chemical potentials 
are estimated as follows. The confinement energy in the 
QW is approximately twice that in the reservoir. Hence 
the confinement energy causes /xid to rise with respect to 
/j,2D ■ The ratio of kinetic energy in the QW to that in the 
reservoir is £fi/£f2 ^ cm. Since the QW is supposed to 
be a ID system, the product an must be small. Hence the 



Fermi energy epi in the QW is noticeably smaller than 
the Fermi energy ep2 in the reservoir. This results in 
lowering /xid relative to /i2D- The ratio of the exchange- 
correlation energies in the case of rs < 1 is estimated as 
£xci/£xc2 ~ (an)~^/^. The exchange-correlation energy 
in the QW is seen to be larger than that in the reservoir. 
Taking into account that the exchange-correlation energy 
is negative, we conclude that it lowers /l(id with respect 
to ^2D- 
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FIG. 1. Dependence of the chemical potential difFerence 
A/x between the decoupled QW and the 2D electron reser- 
voir on the wire radius a (in units of the Bohr radius) for 
various density parameters: (a) rs < 1, (b) > 1. The cal- 
culations were done with the following parameters of the 2D 
layer: d = 2a, Rs = 0.5. 

If the wire is now coupled to the 2D reservoir, we can 
distinguish three cases. For A/i > 0, electrons are trans- 
ferred from the QW to the reservoir, so that the positive 
background charge dominates the electron charge in the 
wire. If A/1 < 0, the reservoir supplies electrons to the 
QW producing an excess negative charge. When A/i = 0, 
the electron dpisity is not redistributed between the wire 
and reservoir .t3 

The electron density redistribution continues until an 
equilibrium state is attained in the whole system such 
that there exists a uniform electrochemical potential. As 
this takes place, a charge and a built-in electric field Eq 
appear in the QW. If the reservoir conductivity is high, 
the built-in field satisfies the condition e /g dx Eo{x) = 
A/i. Such a charge transfer is similar to what happens 
when a contact potential difference appears in classical 
3D systems. However, an essential difFerence is in the dis- 
tance over which the contact field is screened. In the 3D 
case the contact field is screened over a finite length (e.g., 
the Debye length or the Thomas- Fermi screening length). 
In the mesoscopic structure considered here, the contact 
field is produced by the charge, one part of which is sit- 
uated in the ID wire and the other on the surface of the 
reservoir adjacent to the wire. One can easily see, that 
owing to the 3D nature of the electric field, any distribu- 
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tion of charges in a QW cannot screen the contact field 
over a finite distance. Thus the question arises of how 
the charge density and the electric field are distributed 
in the QW and over the reservoir surface. In the case of 
a high enough electron density, the interaction effects are 
not strong, and this problem can be soted analytically 
using the Thomas-Fermi approximation.E3 In the present 
paper, we study this problem using a numerical solution 
of the Schrodinger equation coupled to the Poisson equa- 
tion within the Hartree-Fock approximation. 



III. MODEL 



Finding the self-consistent electron density and the po- 
tential in a QW coupled to electron reservoirs is a rather 
complicated problem because the electron density redis- 
tribution between the QW and the reservoirs produces 
a strong variation of the electron energy in the QW. 
The energy variation is estimated as A^, which is shown 
above to be of the order of the Fermi energy. In essence, 
the QW and the reservoirs should be considered as a 
unique quantum system. To our knowledge, such a prob- 
lem has not been studied to date. In the present paper, 
we investigate it using a simplified model based on the 
Hartree-Fock approximation, which allows one to take 
adequately into account the charge accumulated in the 
QW both for the equilibrium state and under the far 
from equilibrium conditions appearing when an exter- 
nal voltage is applied. This approach allows one also to 
study the exchange interaction effect on the charge accu- 
mulated and the electric potential. However, it does not 
take into account the electron correlation energy. The 
ratio of the correlation energy to the exchange energy 
depends on the electron density. For ID conductors this 
ratio can be estimated according to Ref. 
relation energy is negligible when Ts < 
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The cor- 
1. Under this 
condition the Hartree-Fock approximation is justified. 

It is instructive to estimate numerically the number of 
electrons in a GaAs QW when rg = 0.5. The Fermi en- 
ergy and the electron density are respectively £f ~ 4 meV 
and n « 6 X 10^ cm~^ . Our computation procedure works 
well when the QW length is not too large in comparison 
with the Fermi wavelength. If one puts the QW length 
equal to 0.3 /xm, the total number of electrons in the 
QW is estimated as about 18. Thus the system contains 
about ten electrons in a QW open to reservoirs. The 
number of electrons really existing in the wire and their 
density distribution is determined by the QW length, the 
background charge, the applied voltage, and the e-e in- 
teraction energy. The parameters of the QW estimated 
above are easily realizable in experiment. 

In order to investigate cases with different relative po- 
sitions of the chemical potentials in the QW and reser- 
voirs, we introduce a positive background charge density 
erib in the QW, which is considered as a parameter of our 
model. By varying en;, it is possible to realize any rela- 



tive position of the chemical potentials of the uncoupled 
system. Charges on the QW surface are not taken into 
account in our present consideration. 

Another simplification is that only the lowest subband 
in the QW is considered. 

A ID sketch of the energy diagram of a QW with leads 
and the electron flows (in a far from equilibrium situa- 
tion) is shown in Fig. 0. Here Uq is the confinement 
energy in the QW and fi± are the chemical potentials in 
the reservoirs. The positions of and /i_ relative to the 
conduction band bottom of the corresponding reservoirs 
are fixed because the electron gas in the reservoirs is in- 
compressible. However, and fj,- are shifted relative 
to each other in the presence of an applied voltage Va . 



injected electron flow 




FIG. 2. Energy diagram of a QW with leads and the elec- 
tron flows. 

The following energies contribute to the potential 
shape of the structure: the confinement energy in the 
QW; the potential produced by the external voltage 
source; the Hartree and exchange energies in the QW. 
In the reservoirs > L/2, L being the wire length) 
the electron density is taken to be so high that the e- 
e interaction energy can be ignored in comparison with 
the kinetic energy. Moreover, the reservoirs are assumed 
to be ideally conducting and hence can be treated as 
equipotentials. Inside the QW the e-e interaction poten- 
tial U{x,x') is determined by both the direct interaction 
of electrons with each other and the indirect interaction 
via image charges induced by electrons in the reservoirs. 
Because of this, the interaction potential depends on the 
coordinates cc, x' of the interacting electrons separately, 
rather than simply on their difference. In the QW the 
one-electron wave functions ipr,k^.s{x) are characterized 
by quantum numbers: r = ±l(r = +l indicates elec- 
trons incident on the QW from the left reservoir and 
r — —1 corresponds to electrons incident on the QW 
from the right reservoir); kj. is a wave number in the left 
(r = -fl) or right (r = —1) reservoir; s is a spin vari- 
able. In this paper we ignore changes in s and suppose 
that the states with opposite spins are equally occupied. 
The effects of possible spin polarization in a QW will be 
considered elsewhere. 
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The electron transport in the reservoirs, close to the 
transition between the reservoir and the QW, is modeled 
by a ID Schrodinger equation without interaction. This 
model allows one to simulate adequately the transmis- 
sion probability between the reservoirs and the QW. The 
value of the transmission probability calculated in this 
way as compared with that for a true 2D to ID tran- 
sition at the interface of the reservoir and the wire has 
been investigated by direct comparison ofJiie transmis- 
sion probabilities of these two geometries. EZl The differ- 
ence depends on the energy but does not exceed 15% even 
close to the transmission threshold. The physical reason 
for the success of the purely ID model is that near the 



chemical potential level, of all the electron waves in the 
reservoir only those couple effectively to the wire that 
have a wave vector that is nearly parallel to the wire 
axis (collimation effectEa). We also emphasize that in 
the problem under investigation the electron density of 
the ID Schrodinger problem is needed only within the 
wire. It is unimportant for the reservoirs because they 
are treated as equipotentials and therefore the calculated 
electron density distribution in the reservoirs does not 
directly affect the e-e interaction in the QW. The calcu- 
lated electron wave functions in the reservoirs affect only 
the transmission probability through the contacts. 
In the reservoirs, ^r,K,s{x) is thus 



, , . _ \ eyi-g[ikr{rx + L/2)]+ Rre-xp[—ikr{rx + L/2)] if rx<—L/2 , , 

n'r,KAx) -yxr eMKirx - L/2)] if rx > L/2 , 

where {k'^Y = kf + r x 2meValf^ , and is the applied voltage. 
In the QW, 'ipr,kr,s{x) is determined by the equation 

~2m ^ U(){x) -Uc^t{x) +Uh{x) + V'r,fc„s = £r(fcr)V'r,fe,,s • (5) 

In Eq. (H) the potential energy has the following components. Uq{x) is an effective potential that simulates the 
electron confinement in the QW. In the simplest case, we can assume that Uo{x) = Uq = const for |a;| < L/2 and 
Uo{x) = for \x\ > L/2. Ur is the Hartree energy, 

L/2 

Uh{x)= J dx'U{x,x')[n{x')^nb] , (6) 

-L/2 

with nt being the positively charged background density, n{x) being the electron density. 
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r=±.s-' 



(7) 



and f{kr) the electron distribution function in the reservoirs. In Eq. we assume for simplicity that the radial 
component of the background charge density is the same as the electron density. 

The e-e interaction potential U{x,x') that appears in Eq. depends on the spatial configuration of the leads. 
In what follows the numerical calculations are carried out for the case where the leads are represented as two plates 
perpendicular to the QW. This configuration is comienicnt for further calculations because in this case a relatively 
simple analytical expression is obtained for U{x, x')E3t3 This form of the interaction potential allows one to take into 
account not only the direct Coulomb interaction of electrons but also their interaction via image charges induced on 
the lead surfaces. The interaction potential is 

rj(^ ^2 7 dy 2 r sinh[y(l/2 + 0]sinh[y(l/2-e')] if C < 

eLy sinhyl^^l \ sinh (1/2 - 0] sinh [y (1/2 + ^')] if ^ > ^' , 



where ^ = x/L and Xy is the Fourier transform of the radial density, which is taken to be Xy — 6xp [— (aj//2L)^] . 
Using the analytical expression for U{x,x'), instead of direct solution of the 3D Poisson equation, highly facilitates 
computations. 

Hex is the exchange energy operator, 

L/2 

Hcx{x)lpr,k,.,s = J dx'U{x,x')ncxix,x')lpr,k^^s{x') , 
-L/2 
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where 



CXD 



For the reservoir configuration that we consider here, the 
"external potential" is a linear function of x: Ucxt = 
eVa{x/L + l/2). 

The energy er(fcr) in Eq. (j^) is expressed in terms of 
the wave vector and the applied voltage Va, 



2m 



eVaSr,- 



where we assume that the energy reference is fixed at 

X — —CO. 

In addition, we require continuity of the wave functions 
determined by Eqs. (^) and (|4|) and their derivatives at 
the reservoir- wire interfaces x = ±L/2. The distribution 
functions /(fc^) in the reservoirs are taken in the form 
of the Fermi functions with the temperature T consider- 
ing the fact that the Fermi level in the right reservoir is 
shifted down by eVa with respect to the left one. 

Inside the QW no distribution functions are assigned. 
The electron distribution over the energy is determined 
by the electron flows from the left and right reservoirs 
and the interaction processes inside the QW. The exter- 
nal voltage produces a variation of the electron flows, as 
a consequence of which the electron density is changed. 
Ultimately, this results in the self-consistent variation of 
both the electron states and their occupation for both 
the left and right moving particles. 

The wave functions that we consider in this work 
characterized by a continuous quantum number kr 
Hence, tjj should be considered as a function of two vari- 
ables X and kr- Equation (||) is an integro-differential 
equation with respect to the variable x and an integral 
equation with respect to the variable kr- We develop a 
numerical scheme for the solution of this equation on 
a grid spanning the two variables. The computation 
method is described in the Appendix. It is worth not- 
ing that in the case where the voltage is applied, the 
wave functions are found without using any expansion in 
terms of the undisturbed wave functions. 

The numerical computations were performed using the 
32-processor computer system Parsytec CC. 



IV. THE EQUILIBRIUM STATE 

First, we consider the equilibrium state that appears 
in the absence of an applied voltage Va = 0. In order to 
realize the three cases {Afj, > 0, A/x < 0, and A/i = 0) 
described in Sec. ||, we vary the density of the positive 
background charge ni,. In doing this it is convenient to 
compare ni, with the characteristic density 



where fj,o is the equilibrium level of the chemical potential 
in the system. This quantity has a simple physical mean- 
ing in the case where the exchange and correlation inter- 
action is absent. It is the background density that deter- 
mines which of the three cases is realized in the Hartree 
case. If ni, ~ rip, electrons are not redistributed between 
the QW and the leads in the equilibration process, if 
rib > no, electrons flow from the QW to the reservoirs, 
and if rib < no-, electrons are transferred from the reser- 
voir to the QW. Of course, turning on the exchange and 
correlation interaction shifts the value of the background 
density at which electrons are not redistributed. Never- 
theless, as a reference, the value uq remains convenient. 
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FIG. 3. (a) The electron density distribution and (b) the 
potential energy shape in the QW for nt = 1.5no. The dashed 
lines represent the case without exchange interaction. The 
solid lines are obtained by taking the exchange interaction 
into account. The calculations were done for the parameters: 
cm, 

r=0, and Va=Q. 



a = 5 X 10"'^ cm, L/a = 30, Uo = 20 meV, /xo - ?7o = 4 meV, 



The electron density distribution n{x) for the three 
cases is illustrated in Figs. ||, ^ and ^. Shown here are 
also the graphs of the potential energy U{x) that includes 
the confinement energy and the Hartree energy. 



no 



(8) 



U{x) = Uo + U„ 
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but not the exchange energy. The exchange energy is not 
included in U{x) because it is a functional of rather 
than a direct function of x. It is useful to note that Uh 
essentially coincides with the electric potential. 

Let us consider first the case where the exchange inter- 
action is not taken into account. It is illustrated by the 
dashed lines in Figs. ||, ^, ||. In this case U{x) gives the 
full single-particle potential shape in the QW. If ni, > no 
(this corresponds to > 0), Fig. || shows that the po- 
tential shape lies below the Uq energy and hence the in- 
teraction energy is negative. This means that a positive 
charge is accumulated in the QW. It is responsible for 
the appearance of a potential well. 



due to the charging of the QW. 
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FIG. 4. The same as in Fig. ^ but for nt = 0.5no. 

The case A/x < is realized when Ub < no. The 
electron density distribution and the potential shape are 
shown in Fig. ^. The interaction energy is seen to be 
positive and the potential shape in the QW lies above 
Uo- This means that a negative charge is accumulated 
in the QW. It produces a potential barrier that hinders 
electrons in passing through the QW. 

If rift = no , the electron density is essentially not redis- 
tributed between the QW and the reservoirs. However, 
Friedel oscillations of the electron density appear near 
the contacts, Fig. ||. The Friedel oscillations are also 
observed if there is carrier transfer, but they are super- 
imposed on the much stronger variation of the potential 
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FIG. 5. The same as in Fig. | but for Ub = l.Ono. 

Let us now consider the exchange interaction ef- 
fect. The electron density distribution and the potential 
shape, calculated by taking into account the exchange 
interaction, are shown in Figs. |[ ^ || as solid lines. The 
exchange interaction is seen to result in an increase of the 
electron density. This is a consequence of the fact men- 
tioned in Sec. || that the exchange interaction decreases 
the chemical potential in a QW. That is why more elec- 
trons come into the QW when the equilibrium state is 
established. Correspondingly, the negative charge in the 
QW increases, which results in the growth of the energy 
U{x). However, this does not yet mean that the elec- 
tron states with energy lower than U{x) are necessarily 
states decaying in the QW, since the exchange interac- 
tion lowers the effective barrier between the QW and the 
reservoir. 

Since it would be incorrect to consider the exchange 
energy as a function of x, we calculate an average value 
of the exchange energy per particle incident on the QW 
with the energy e(A:-)_) from the left reservoir, 



£^ox — 



(here (• • • ) denotes averaging over the QW length). The 
average exchange energy i?ex is shown in Fig. ^ as a func- 
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tion of e(fc+). The average Hartree energy Eh is also 
given in this figure. It is seen that E^x exceeds Eh for all 
energies. Hence the joint effect of the exchange interac- 
tion and the Hartree interaction consists in an effective 
lowering of the barrier, so that electrons with energy be- 
low the confinement energy Uq can transit through the 
QW without decay. 
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FIG. 6. Average exchange and Hartree energies of electrons 
incident on the QW with the energy e. The calculations are 
for the parameters: a — 5x 10^^ cm, L/a — 30, Uo — 20 meV, 
MO - f^o = 4 meV, Ub = l.Ono, T= 0, and 14= 0. 

Another effect produced by the exchange interaction 
is also seen from Figs. ^, ^, The exchange interaction 
strongly enhances the Friedel oscillation amplitude. This 
result agrees qualitatively with the analytical calculation 
of the interaction effeC|t-|On the transmission through a 
barrier in ID systems.E2l In our case the Friedel oscil- 
lations are generated at the contacts of the QW with 
the reservoirs. Their amplitude in the potential energy 
is quite pronounced but smaller than the Fermi energy 
A*o — Uq in the QW. For the discussion which follows it 
is important to remark that the Friedel oscillations are 
superimposed on the smooth variation of the potential 
produced by the charge accumulated in the QW. Even 
if the exchange interaction is fully taken into account, 
this smooth component has an amplitude that is larger 
than that of the Friedel oscillations. Due to the smooth 
variation of the potential the QW becomes nonuniform. 



V. LINEAR CONDUCTANCE 

The model that we have developed above allows us to 
find the electric current arising when an external volt- 
age is applied. The current is calculated as the sum 
of the partial currents of the states ^r,kr.s, taking into 
account their occupation. The linear regime is realized 
when eVa <SC (/^o ~ Uq). In this case we have obtained the 
dc conductance as a function of the chemical potential 
/io- The results of these calculations for zero tempera- 



ture and for several densities of the background charge 
rib are given in Fig. 0. 




A£p = \i- f/y (meV) 

FIG. 7. Linear conductance as a function of the chemical 
potential ^o- Lines 1, 2, and 3 present the calculations within 
the Hartree approximation for several background densities: 
nb/no = 1.5, 0.5, 1.0, respectively. Line 4 is obtained by in- 
cluding the exchange interaction for Ub = uq. The parameters 
used in the calculation are a = 5 x 10~^ cm, L/a — 20, 
10 meV, r= 0, and Va= 0.01 mV. 

The conductance oscillations with varying chemical po- 
tential are a consequence of the nonadiabatic reservoir- 
wire interface. The rapid variation of the potential at this 
interface leads to backscattering and, if the electron wave 
is coherent over the entire wire length, to resonances. 
The oscillations have the same origin as the resonances 
observed in over the barrier transmission of noninteract- 
ing particles. A similar effect also appears in the trans- 
port of noninteracting electrons through a narrow, .-ballis- 
tic, nonadiabatic constriction in a 2D electron gas.ES Our 
calculations show that the e-e interaction changes the ef- 
fective potential barrier that electrons have to overcome 
in passing from one reservoir to another. 

First, we discuss the results obtained within the 
Hartree approximation, when exchange interaction is ne- 
glected. These results are represented by the lines 1-3 
in Fig. for several different background densities m- 
With increasing rit, the oscillations become more fre- 
quent, which means that the effective wave number of 
the electrons is increased. Exactly the same behavior is 
demonstrated in Figs. |[ ^ |^, namely, with increasing rib 
the potential in the QW is shifted downward causing the 
kinetic energy to increase. 

A similar effect occurs when the exchange interaction 
is turned on. It is demonstrated in Fig. [7| by curve 3 
(obtained by ignoring the exchange interaction) and by 
line 4 (obtained by including exchange interaction), the 
background charge being the same in both cases. The 
exchange interaction is seen to make the conductance os- 
cillations more frequent. The reason for this effect is that 
the exchange interaction results in an effective lowering of 
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the potential energy of the electrons and consequently in 
an increase of their kinetic energy. In order to assess the 
exchange interaction effect on the effective potential, it 
is instructive to see how the exchange interaction, affects 
the spectral density of electrons, i. e., IV'fcP integrated 
over the QW length. This is illustrated in Fig. ||. The ex- 
change interaction allows the electrons with energy below 
the confinement energy Uq to pass through the QW. 
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FIG. 8. Exchange interaction effect on the spectrum of the 
electron density {4>k\tpk}- Line 1 presents the spectral density 
obtained by including the exchange interaction; line 2 is the 
spectral density without the exchange interaction. The cal- 
culations are for: a = 5 x 10~^ cm, L/a — 30, Uo= 20 meV, 
MO - [^0 = 4 meV, Ub = l.Ono, T= 0, and 14= 0. 

It is interesting to note that, despite the fact that 
Friedel oscillations are present in the QW, no noticeable 
suppression of the conductance is observed. The effect of 
conductance suppression by a periodic potential associ- 
ated with Friedel oscillations was considered for infinite 
ID systems with a i5 potential in Refs ^ and This 
phenomenon is connected with the fact that a periodic 
component of the potential suppresses the transmission 
of the electrons with energy near the Fermi level across 
the QW (a gap appears at the Fermi level). The absence 
of this effect in our system is a consequence of two facts. 
First, the QW has a finite length. Second (and no less 
essential) , the Q W becomes inhomogeneous owing to the 
electron density redistribution between the QW and the 
leads. As a consequence the kinetic energy at the chem- 
ical potential level and the Friedel oscillation period be- 
come dependent on the position in the QW. This is why 
the resonant interaction of electrons at the chemical po- 
tential level with the Friedel oscillations is destroyed and 
the electron passage is not suppressed. 



VI. NONLINEAR TRANSPORT AND 
MULTISTABILITY 

A significant redistribution of the electron density be- 
tween the QW and the reservoirs occurs under far from 
equilibrium conditions when the applied voltage exceeds 
the Fermi energy. Electrons are injected from the left 
reservoir (cathode) while the electrons entering the QW 
from the positive reservoir (anode) are scattered back 
inside the QW. As a consequence, the electron density 
decreases in the QW (roughly speaking to one-half of 
the equilibrium density) though the positive background 
charge is unchanged. Because the positive charge is dom- 
inant, a potential well appears in the QW, with the po- 
tential shape being distorted by the external potential, as 
illustrated in Fig. ||. Therefore the kinetic and potential 
energies are greatly changed. The change in the potential 
energy produces variations in the wave functions (includ- 
ing even a possibility for resonant states to appear) and 
the electron density distribution. In this way feedback 
arises between the electron density and the potential in 
the QW, which is an important mechanism in nonlinear 
transport. It is that mechanism which is realized in the 
model proposed. 

A complete numerical analysis of the nonlinear trans- 
port properties in the wide range of applied voltages 
within the Hartree-Fock approximation meets some dif- 
ficulties caused by the long computation time. In this 
paper we restrict our consideration to the Hartree ap- 
proximation, which is reasonable at high enough voltage 
because the exchange energy decreases when the kinetic 
energy of electrons is increased^ The calculations were 
carried out using the method of pseudotime evolution to 
the steady solutioiO described in the Appendix. It turns 
out that in some range of applied voltage an instability 
of the evolution process appears. The instability origin 
is not connected with the computation process but is 
caused by real behavior of the system. 

The mechanism of the instability is as follows. When 
the applied voltage is high enough (compared to the 
Fermi energy), the electron flow injected from the nega- 
tively charged reservoir is the only flow in the QW. Let a 
velocity fluctuation appear in some portion of the wire. 
To be definite, let us assume that the velocity is increased 
above its stationary value. Since the total electron flow is 
limited by the contact, it is not disturbed by this fluctu- 
ation. Hence, the continuity of the current requires that 
the electron density decreases. This leads to a growth 
of the positive (net) charge, because electrons cannot 
completely neutralize the background charge. The ex- 
cess positive charge causes the potential energy of the 
electrons to decrease. Under the condition of ballistic 
transport, this results in a new increase of the veloc- 
ity, and so on until some nonlinear process stabilizes this 
instability. In our model this is achieved by a redistri- 
bution of the overall electron density and a reshaping of 
the potential distribution in the QW. In such a way the 
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potential shape is switched from one state to the other 
under the condition that both states are characterized by 
the same potential difference across the QW ends. 
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FIG. 9. Dependence of the average kinetic energy of elec- 
trons on the apphed voltage. The calculations are for the fol- 
lowing parameters: a = 5 x 10~^ cm, L/a = 20, Uq— 10 meV, 
Mo — t^o = 3 meV, rib/no— 1.0, and T= 0. 

In order to describe the transition from one shape to 
the other as a continuous process it is necessarily to char- 
acterize the nonequilibrium state of the system by a pa- 
rameter other than the applied voltage. This parameter 
should distinguish states with different potential shapes 
and the same potential difference across the ends. As 
such a parameter, the mean kinetic energy £^kin of elec- 
trons in the QW can be used, 

_ Y.r=±S dkr{kr\f\kr) 
-^kiii — i2 •) 

where T is the kinetic energy operator. 

This conclusion is similar to what is known in the the- 
ory of hot electron instabilities in semiconductors. The 
electron heating by the electric field results in S- and 
N-shaped current-voltage characteristics. Under these 
conditions, the electron temperature uniquely determines 
the state of the system.E3 In our case the kinetic energy 
of the electrons is a direct analog of the electron temper- 
ature. 

We have developed the algorithm that allows one to 
solve our problem in the case where the mean kinetic en- 
ergy is fixed rather than the applied voltage. This algo- 
rithm is described in the Appendix. In this computation 
scheme all quantities (including Va) are determined by 
^'kin- We have found that this algorithm gives stable re- 
sults and the nonequilibrium state of the system under 
investigation is uniquely determined by -Ekin- In partic- 
ular, the dependence of E'kin on the applied voltage is 
shown in Fig. ^. The kinetic energy is seen to have sev- 
eral values for a given voltage Va, while Va is uniquely 
defined by i^kin ■ Correspondingly, several nonequilibrium 



states, with different distributions of the electron density 
and the potential, are possible at a given voltage. It is 
obvious that the states are not all stable with respect to 
time-dependent fluctuations. 

The multistability, and in particular bistability, phe- 
nomenon described above may be useful for understand- 
ing the negatiiig differential conductance observed in 
quantum wires.l23 



VII. CONCLUSIONS 

In this paper a QW of a finite length coupled to reser- 
voirs is considered as a unified system. We have found 
that the electron density is substantially redistributed 
between the wire and the electron reservoirs when this 
system is formed. As a consequence of this process, a 
QW can acquire a net charge. The charging of the wire 
is caused by the chemical potential difference between the 
wire and the leads which exists if these two subsystems 
are decoupled. This phenomenon is similar to the contact 
potential difference in a classical multiconductor system. 
The structures of the charge density and the electric po- 
tential distributions differ according to the chemical po- 
tential difference: (i) a positive charge is accumulated in 
the wire and a potential well is developed there; (ii) the 
wire is charged negatively giving rise to a potential bar- 
rier; (iii) the wire remains uncharged as a whole. In all 
cases, Friedel oscillations are present which are generated 
at the nonadiabatic contacts of the QW with reservoirs. 
They are superimposed on the relatively smooth profile 
of the potential produced by the charge accumulated in 
the wire. This smooth potential has a large amplitude. 
The Friedel oscillation amplitude is strongly enhanced if 
the exchange interaction is included. 

Variation of the electron density in a QW due to elec- 
tron redistribution between the wire and the reservoirs 
produces a significant effect on dc conductance. This ef- 
fect is connected with the change of the kinetic energy 
of the electrons due to two factors: (i) the variation of 
the accumulated charge and the potential variation asso- 
ciated with this charge; (ii) the variation of the exchange 
energy. The exchange energy substantially lowers the ef- 
fective potential barrier that electrons have to overcome 
when passing from one reservoir to another. 

The Friedel oscillation potential does not suppress the 
conductance because the QW becomes inhomogeneous as 
a consequence of the electron density redistribution be- 
tween the QW and the leads. The electron density vari- 
ation due to the voltage applied across the leads and the 
change of the effective potential shape, associated with 
this variation, in the wire is an important mechanism for 
nonlinear transport. 

The most interesting consequence of the charging ef- 
fect in the wire is the instability that arises under a high 
enough applied voltage. The instability shows itself as a 
spontaneous increase of the kinetic energy of the injected 
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electrons at a given applied voltage. In a QW with a 
nonadiabatic connection to the reservoirs, the develop- 
ment of instabilities results in the appearance of multi- 
stable states, i.e., in the existence of several stable states 
at a given voltage. 

We conclude by emphasizing that the charging effect 
analyzed here is a general phenomenon that might be im- 
portant in many mesoscopic systems containing electron- 
ically different compoTOds. Examples of current interest 
are carbon nanotubes,EZI hybrid normal-superconducting 
systems, and atomic quantum point contacts. 
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APPENDIX A: THE COMPUTATION METHOD 

The problem of finding the wave functions for a QW 
coupled to reservoirs can be reduced to the solution of 
Eq. (^ in the inner region — L/2 < x < L/2 and match- 
ing ipr,k,s and its derivative at the boundaries x = ±L/2 
to the wave functions (^) in the outer regions. Com- 
bining the matching equations for ipr,k,s and ip'r k one 
can exclude the coefficients Rr and Tr to get finally the 
following nonuniform boundary conditions for ilJr,k,s- 



rx = -L/2 
rx = L/2, 



(A.l) 



where kr and fc(, are defined in Eq. (Q) 

The wave functions ipr.k^.s are considered as functions 
of two continuous variables: the space coordinate x and 
the energy = h^k^/2m. Thus ipr,kr,s = '^r,s{x,er), 
where x and Sr are varied respectively in the regions 
—L/2 < X < L/2 and < < £m, with the upper 
boundary em being equal to /ig -I- iksT {ks is the Bolz- 
mann constant, T is the temperature). In this region the 
uniform grid {xi,ej}i=Q^Ni;j=Q,N2 is constructed. 

The equation for ipr.k.sixi, ej) on the grid is obtained 
from Eq. (|^) with use of the integro-intcrpolative method 
for the node presentation on Xi and the trapezium for- 
mula when calculating the integrals on kr for the nodes 



on Ej. This results in a set of nonlinear finite-difference 
equations that can be symbolically presented in the form 



M[^]^ = F , 



(A.2) 



where is the wave function vector to be found and 
M[^] is the nonlinear operator. The matrix equa- 
tion ( A.2 ) is inhomogeneous as a consequence of the 
boundary condition s (lA.lh . 

The equations ( |A.2| ) are solved by the iteration 
method. However, the commonly used successive approx- 
imation scheme (M[*(')]^'('+i) = F, with I being the it- 
eration number) turns out to be badly convergent. We 
use thoinethod of pseudotime evolution to the steady so- 
lutionEJ More specifically, we use the two-layer iteration 
scheme of this method. In this scheme the approximat- 
ing matrix M^'^ is introduced, which is calculated with 
the use of the iteration process 



n 



M[^'(')] - M«, I = 0,1,2, 



(A.3) 



where r; is a pseudotime parameter. The choice of r/ al- 
lows one to attain the best convergence of the iteration 
process. As the starting value of M'^°) we use M[\E' — 0], 
i.e., the M matrix for noninteracting electrons. During 
the iteration process, ^^'^ is calculated with the use of 
the equation 

]V[(Ovj/(0 ^ F ^ 

The pseudotime r/ is determined by the M*^'^ operator 
spectrum. The optimal convergence is attained when 



Amin(M(0)+A„,ax(M(0) ' 



where Amin and Amax are the lowest and highest eigen- 
values of M*^'^. The iteration process is ended when the 
following condition is fulfilled 



max 



,(0 



< <5 



where to,-'-* is an element of the M^'^ matrix. In the 



present paper, 6 was chosen to be 10~^. 

The above method is successful when the system un- 
der investigation has a unique solution. However, at some 
fixed values of the applied voltage the computation shows 
an instability. In the course of the pseudotime evolution 
process the calculated quantities (such as the potential, 
the electron density, the kinetic energy) are randomly 
switched between several values. This is connected with 
the fact that the state of the system is not uniquely de- 
termined by the calculation scheme where the applied 
voltage is fixed. 

A unique description of the system is achieved by us- 
ing the mean kinetic energy i?kin of the electrons as the 
parameter that defines the nonequilibrium state of the 
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system. We have developed computation algorithm that 
allows one to vary iJkin continuously, in other words, we 
solve the problem using E'kin as the fixed parameter in- 
stead of the applied voltage. 

An essential question appearing in this algorithm is 
how the applied voltage 14 should be defined when Skin is 
given. The equation defining Va is obtained from Eq. (|^). 
Multiplying this equation by ip^ integrating it over 
kr and over x, and summing over r one gets an equation 
of the following form: 



(A.4) 



where A [4"] and S[\E'] are functional of the electron wave 
functions. Solving this equation with respect to Va one 
gets Va as a functional of with Skin being a parameter. 



Va = <^E,J^] 



(A.5) 



W hen solving the problem with Skin as a parameter, 
Eq. (|A.5D should be taken into account together with 
Eq. ( |A.3D . This system of equations is solved using the 
above pseudotime evolution method and two-layer iter- 
ation scheme. The set vi'"* approximating Va is defined 
as 



Tl 



= -yj') , ; = 0,1,2,... . 



As the starting value of the yi''' set, the arbitrary value of 
Va in the stability region close to the instability threshold 
can be used. In this generalized procedure the pseudo- 
time Tl is chosen taking into account the spectral prop- 
erties of the total matrix M^') ® vi'' . 
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